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Abstract 

Perturbations in a universe with matter and cosmological constant are computed in 
terms of the primordial scalar perturbation up to second order. The calculation is easy 
- a generalization of the spherical dust solution. Our expression might allow accurate 
measurements of the primordial non-Gaussianity. 
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I. INTRODUCTION 

Inflationary models ^ explain the origin of cosmological perturbations, and predict the 
primordial super-horizon non-Gaussianity which is small in the standard one-field slow-roll 
case 12]. A simple and natural extension of the standard model, however, predicts a much 

n n 

larger non-Gaussianity [3|, 1^- To measure the primordial non-Gaussianity one needs to 
compute the current sub-horizon perturbations in terms of primordial super-horizon pertur- 
bations up to second order in the amplitude. This is done here for the universe filled with 
matter and cosmological constant. The pure matter case has been computed by in that 
paper we also discuss observations of non-Gaussian effects in galaxy surveys. 

Neglecting tensor (gravitational) perturbations, we assume that primordial perturbations 
are pure scalar. Following [3, ^B], we use a gauge invariant variable ( that remains constant 
outside the horizon to characterize the primordial perturbations. On super-horizon scales ( 
is defined as follows: is proportional to the local scale factor measured on uniform local 
Hubble parameter hypersurfaces. 

The degree of non-Gaussianity of primordial perturbations /„; can be defined by assuming 
that ( is of the form ^ 

C = g- {3/5)fnig', (1) 

where g is Gaussian. The standard one-field slow-roll inflation predicts /„; < 0.05 
Other inflationary models predict |/„/| > few. The WMAP results are —58 < fni < 134 

Sky Survey should provide similar accuracy, and 
11|. 

The non-Gaussianity of sub-horizon density perturbations is due to Newtonian evolution, 
relativistic effects, and primordial non-Gaussianity. As we will see, the ratio of relativistic to 
primordial contributions is ~ l//n«; accurate determination of fni must include relativistic 
effects. 

In §2 we derive the second-order expression for the density perturbation in terms of (. 
In §3 we generalize the standard Newtonian calculation of the density skewness. 



at 95% confidence 'lO'|. The Sloan Digital 
future galaxy surves can be much better 
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II. SECOND-ORDER DENSITY PERTURBATION IN TERMS OF C- 



Results 



Consider flat matter plus cosmological constant universe, with no primordial gravitational 
waves. We use the comoving synchronous gauge. The first-order metric is 

ds"^ = dt^ - ao(t)(l + 2C)(5»j + 2rCij)dxidxj. (2) 

Here ao(t) is the unperturbed scale factor; C,{x) is the primordial perturbation, which is a 
function of x only; indicates spatial derivatives; T{t) is the linear growth factor given in 
the Appendix. For zero cosmological constant, A = 0, we have r oc t^/^ - the familiar linear 
growth of perturbations in the matter dominated universe. 
The density perturbation to second order is 

5 = -rCu + y ((1 + + (1 - C)CS) + I (4Ca - C^) . (3) 

Here C{t) is a dimensionless numerical factor given in the Appendix. For A = 0, C = 3/7 
and we recover j^. 

The first two terms in Q represent the second-order Newtonian perturbation in the 
synchronous coordinate system 



121 |l3[ . The third term represents the relativistic effects. 
The first-order metric can be derived along standard lines jlj]. The second order density 
perturbation is obtained in the next subsection. The applications of this expression are 
discussed in §3 and §4 . 

B. Derivation 

We calculate the second order density perturbation under spherical symmetry. Then we 
simply generalize to the non-spherical case. 

Consider the well known exact solution of the Einstein equations - the spherically sym- 
metrical dust with cosmological constant. The metric is 

ds^ = de - -^^^^dx^ - a'x'ide^ + sin^ Odcj)^). (4) 
1 — x^ I 



The matter density is 

^ 87TG{x^a^y 

Here R{x) and F{x) are arbitrary functions of x but not of t, while a{x,t) and p{x,t) are 
functions of both x and t, and a satisfies the Friedmann equation 

aa'^ + a/R^ = F + Aa^S. (6) 

Dot and prime are time and x derivatives. 
We can choose 

F oc e^^, (7) 
then at t — equation © gives a oc e^t^/^. If we also choose 

= --C - C'\ (8) 

X 

equation (@)) gives the early time (superhorizon) metric of the form 

ds^ = dt^ - al{t)e^^dx^, (9) 
where anft) oc t^/^, t — > 0, is the unperturbed scale factor. We see that what we have called 

nn 

( is indeed the standard gauge invariant quantity that inflation theories like [2, |3|| claim to 
predict to second or arbitrary order. 

Having expressed the two free parameters of the spherical model, R{x) and F{x), in terms 
of C{x), we can now calculate the spherical density perturbations in terms of ( to second 
order. The dimensionless density perturbation 6 is defined by 

p={l + S)po, (10) 

Po is the unperturbed matter density. It is convenient to define 6a by 

a = e'^aoit){l + 5a). (11) 

Then equation © gives 

6 = -36a - x6'a + 66l + ixSaS'a + xX + ^'CK, (12) 
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and equation (jH)) gives 

5a = 6al{l-C6al-2C+]^xC'). (13) 

where 

5al = T-C'. (14) 
X 

Finally we put into p2|l and get the spherically symmetrical density perturbation 
to second order 

5 = -r(C" + -CO + r^C" + ^±^CC" + + ^(4CC" + -CC - C'^). (15) 

X X X^ I X 

To get the corresponding expression without the spherical symmetry assumption, we 
write down the most general expression for a scalar which contains the same powers of C, r 
and V as the expression (fT3j): 

b = -tC,Q, + T^C^Cl + C,Cl) + r{C,CC^^ + C,C^)- (16) 

Now assume that C, is spherically symmetrical. Then ()16|) should reproduce ()15|) . We have 
5 arbitrary coefficients in (fT^ in terms of which we should get the 8 known coefficients of 
(fTK|l . If we could do it, it would be a check of our spherically symmetrical expression. We 
can, and the result is given by (jHI). 



III. THE SKEWNESS OF THE DENSITY PERTURBATION 

The skewness of the density perturbation S =< 5^ > / < 5^ >^. The computation should 
include the Jacobian of (2), J = 1 + 3C + tQi. To lowest order, 

S - (4 + 20 - -L|(;„, , |) (n, 

Here A;^ =< > / < (f >. For A = 0, neglecting primordial and relativistic contributions 
of the second term, we obtain the Newtonian result 4 + 2C = 34/7. 

The actual perturbation spectrum is scale-free, and our formula is not directly applicable. 
However it may serve to compare primordial and relativistic contributions to the sub-horizon 
non-Gaussianity. 
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IV. CONCLUSIONS 



We have expressed the large-scale cosmological perturbations in terms of the primordial 
scalar perturbations up to second order. Our expression may allow to measure primordial 



non-Gaussianities using galaxy surveys 

I thank Raisa Karasik and Roman Scoccimarro for useful discussions, and P.J.E. Peebles 
for telling me about the spherical perturbation case. This work was supported by the David 
and Lucile Packard Foundation. 



APPENDIX A: r(t) AND C{t) 

The linear growth rate in the synchronous comoving gauge coincides with the Newtonian 
expression for perturbations in a universe with matter and cosmological constant. We just 
need to re-derive it by perturbing eq. (jH)) according to (jllj) . To first order we get 

1 C' 

aoSai + AnGpoaoSai = . (Al) 

ao X 

We write the solution in the form (jl4j) . Then 

Tit) = Dit) / , , (A2) 

where D{t) is another growth factor (the homogeneous solution of ()A1|) ): 

\nD{t) = -A.G r ^^1/^0(^1)^0(^1) . (A3) 



These equation coincide, after simple rearrangements, with the expressions given in 
On the other hand, the cosmological perturbation theory ^| gives the growth factor of 
perturbations in the comoving synchronous gauge 

f = a^^ j dtao, (A4) 

which is equal to ()A2I ) . 

At early time the universe is matter dominated, po = l/(67rG't^) and ao oc t'^^^. Then 
eq. ()A3|) gives D = t^^, and eq. ()A2|) gives the familiar growth factor r = ?7^/10 oc t^^^, where 
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7] is the conformal time. At a later time, when the cosmological constant accelerates the 
expansion of the universe, the growth slows down. 

When equation (jU)) is solved to second order, a non-linear growth factor will appear. The 
second-order solution can be written in the form (jl31 114j) . with a dimensionless coefficient 

rm - ^(^) r 2r + 0^x72 + 2aodoDTyD - a^d^rf 

f-l^J — ~^T77\ I "''-i 7\ '■ ' l^^J 

T\t) Jq Daoao 

with all quantities under the integral being at time ti. At early time, C = 3/7 and we recover 
the zero cosmological constant result 0|. 
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